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Abstract. We study the energy properties of a particle in one dimensional semi- harmonic 
rectangular wells and barriers. The integration of energies is obtained by solving a simple 
1/^ , transcendental equation. Scattering states arc shown to include cases in which the impinging 

particle is 'captured' by the semi-harmonic rectangular potentials. The 'time of capture' is 
'(—( ', connected with the dwell time of the scattered wave. Using the particle absorption method, 

^ i' it is shown that the dwell time coincides with the phase time tw of Eiscnbud and Wigner, 

calculated as the energy derivative of the reflected wave phase shift. Analytical expressions 
arc derived for the phase time tw of the semi-harmonic delta well and barrier potentials. 



Key words: exactly solvable potentials; scattering process; resonances; Eisenbud-Wigner 
phase time; dwell time 



^ ' 2010 Mathematics Subject Classification: 35Q40; 35B34; 81U30; 81Q60 

> ■ 

in 

^ . 

^ ■ 1 Introduction 

I One-dimensional models of quantum mechanics are useful in a number of applications in con- 

temporary physics [IllIliaiiiaEllZllSllglllOllIIllIllliaillllSllinil^ often employed 
into the approximations which make tractable the more elaborated three-dimensional systems, 
they allow to get a deeper insight on the physics involved. Their simplicity has made them 
^ ' valuable as academic and research tools. For instance, the concept of effective mass, success- 

^ . fully applied in describing the formation of shallow energy levels due to impurities in crystals, 

leads to one-dimensional systems jl]. Localized lattice deformations can be modelled as a rect- 
angular potential to study isolated transitions, observed in semiconductors, from a bound state 
within a quantum well to a bound state at an energy greater than the barrier height [2]. The 
analysis of resonances O IH El El d E] is transparent for rectangular potentials in either, the 
presence of a background interaction [U [10] , or in free space [HI [121 [El [Ml [El [El [E] ■ Simple 
models of point-like [IHl [l9] , as well as regularized singular interactions [20] , can be obtained 
as limit cases of rectangular potentials |211 \T7\ [22]. The one-dimensional models are also useful 
in the study of supersymmetric quantum mechanics [231 [Ml [251 [211 127] . Based on the Darboux 
transformations [28] (see also j23j ) , the supersymmetric (intertwining or factorization) formalism 
allows the construction of new exactly solvable potentials even if complex energies are involved 
[291 [30l [3T1 [321 [33] . This last property has been implemented to get one-dimensional complex 
potentials behaving as an optical device which both refracts and absorbs light waves [151 133]. 



*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Benasque, Spain). The full collection is available at 
|http://www.emis.de/journals/SIGMA/SUSYQM2010.html| 



2 



N. Fernandez-Garcia and O. Rosas-Ortiz 



Collisions, on the other hand, are modelled as interactions localized in time and space. This 
implies that the involved potential vanishes rapidly enough in space, so that incoming and out- 
going asymptotic states can be represented by wave packets in free motion. The problem is 
usually reduced to the analysis of a one-dimensional effective potential. Main information is 
then obtained from the transmission and reflection amplitudes. Of particular interest, the reso- 
nance phenomenon is experimentally studied in atomic, nuclear and particle physics [H [U [Sj [6] . 
A resonance can be understood as an special result of the scattering process in which the inci- 
dent wave is 'captured' by the scatterer for a while. A measurable time delay (relative to free 
motion) is then associated to the scattered wave [SH [35] . 

The concept of time delay corresponds to the time spent by an scattered particle in the scat- 
tering zone when compared to a free particle subject to the same initial conditions. According 
to Eisenbud and Wigner, this is associated to the energy derivative of the phase shift for binary 
collisions (see [36l EZ])- The existence of a global time delay (as connected to sojourn times), 
and its identity with the Eisenbud-Wigner phase time, have been proved for local potentials 
in [38]. The proof has been extended to an abstract formalism where time delay is defined in 
terms of the expectation values of non- negative normalized functions of compact support [39]. 
In general, the retardation of the scattered wave involves transient effects which are relevant in 
nuclear reactions [401 [41] , and is also fundamental in the characterization of resonances |34[ [35] . 
The time delay is also connected with the Levinson theorem |42j . the density of states in meso- 
scopic conductors [43] and with the study of the photodetachment rate due to weak time-periodic 
electric fields [H]. In one-dimensional systems, the time delay has been studied for step-like 
potentials exhibiting a two-channel structure [45]. Quite recently, it has been shown that time 
delay can be obtained from the eigenvalues of a non-Hermitian Hamiltonian in two and three 
dimensions [l6]. Extensive, clear reviews with complete bibliography can be found in [47[ I48 [ [^. 

The present work is addressed to the analysis of the energy properties of a particle subject 
to the action of a semi-harmonic rectangular potential. This last is either a rectangular well or 
barrier in a background integrated by a free-particle interaction to the right and an oscillator- 
like interaction to the left of the rectangular potential. The model corresponds to a system (the 
rectangular potential) embedded in an environment (the semi-harmonic background), and the 
issue is the study of the modifications on the energy spectrum and resonances of the system 
which are induced by the environment. Thus, the semi-harmonic square potentials are viewed as 
open one-dimensional quantum systems [5]. The bound states and resonances of square poten- 
tials have been studied by using diverse approaches (see, e.g., [HI [151 [16] a-^id references quoted 
therein). In particular, if closer resonances imply narrower widths (i.e., in the case of isolated 
resonances), the transmission amplitude T can be written as a superposition of Fock-Breit- 
Wigner distributions [15] . Then, the position E and width F of a resonance e = E — are 
in correspondence with one of the bell-shaped peaks of T. This result is extended to the case 
of resonances which are not isolated by identifying the position of the peaks with the absolute 
value of e, rather than using Re(e) = E [16]. In the limit where the square potentials become 
a point-like interaction (i.e., a delta barrier or well), the transparency properties involved pro- 
duce no resonances [18]. The situation is different for dipole-like interactions represented by 
derivatives of the delta distribution with respect to the position [17j . 

As open systems, rectangular potentials have been previously studied in a static field as the 
environment [9l[l0]; the reports include the point-like limit [44]. Unlike the rectangular potentials 
in the free-particle background, these last open systems are not isotropic since the environment 
is represented by a potential diverging as itx at = +oo. Thus, though the presence of 
resonances is ensured by adding such a background, the conventional definitions of dwell time 
and time delay are not automatically applicable in these cases. Indeed, most of the approaches 
on the matter consider a potential which is very localized in a finite region of space [47[ 1^8] [l9] 
(see also [38] )• In contraposition, the open systems presented in this paper become zero when 
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they are evaluated at positions to the right of the square potential. This condition represents 
an advantage over the square potentials in a constant field environment since reflection times 
can be now calculated. In this context, the mean time spent by the particle in the interaction 
zone is defined as the dwell time (or sojourn time) involved. However, it must be emphasized 
that the parabolic part of the environment makes not simple the calculation of transmission 
times. Thereby, there is not a clear definition of time delay in terms of dwell time differences 
in arbitrary open balls of the one-dimensional domain of these systems. Instead, the time delay 
can be defined as the difference between the reflection times above mentioned. Here, we shall 
show that the dwell time coincides with the phase time introduced by Eisenbud and Wigner. 

The organization of the paper is as follows. In Section [2] we establish the general expressions 
concerning the semi-harmonic rectangular potentials. The condition for the trapping of particles 
is identified as the finding of zeros of the Jost function. A very simple transcendental equation 
is derived for the numerical integration of energies. In Section [31 the family of semi-harmonic 
rectangular wells of unit area is shown to converge to a semi-harmonic delta well in the sense 
of distribution theory. The corresponding transcendental equation is shown to admit a unique 
root associated to bound states, and the result is compared with the single bound energy of 
the delta well in a free particle background. The resonances obey a rule of distribution in 
the complex plane of the energies that resembles the distribution of the odd bound energies 
of a harmonic oscillator. Section |4] is devoted to the application of the previous results to the 
case of a semi-harmonic rectangular barrier. Interestingly, the limit case of a semi-harmonic 
delta barrier obeys a distribution of resonances which is in correspondence with the even energy 
eigenvalues of the harmonic oscillator. A brief discussion on the Darboux transformations of 
rectangular potentials in a semi-harmonic background is also given. In Section [5] we investigate 
the dwell time of the scattering process associated to the semi-harmonic rectangular potentials. 
This is defined as the time spent by the incident particle in the interaction zone. We use the 
particle absorption method [501 ISTj [52] to determine the characteristic time involved. Then 
we assume that the probability of finding the particle decays exponentially from the moment 
that it is 'captured' by the potential. The dwell time so calculated coincides with the phase 
time Tw of Eisenbud and Wigner. For delta wells and barriers explicit analytical expressions 
of Tw are given. In each case, the time delay tw has local maxima centered at the real part of 
the involved resonances. This last, combined with the results of Sections [3] and [J] for delta-like 
potentials, indicates that the semi-harmonic background induces the real part of the resonances 
to be distributed in correspondence with the energy eigenvalues of the harmonic oscillator. Some 
concluding remarks are given in Section [6) Finally, a short appendix includes some derivations in 
terms of the confluent hypergeometric functions which, although important, can be postponed 
to a later reading. 

2 Semi-harmonic rectangular well 

Consider a particle of energy E which is under the influence of the one-dimensional potential 
x^, X < —a (Region I), 



V{x; a) = < 



-Vo, \x\ < a (Region II), a > 0, Vb > 0. (2.1) 
0, a < x (Region III), 



Potential (j2.ip is depicted in Fig. [TJ The corresponding stationary dimensionless Schrodinger 
equation 

{H - E)^{x) =0, + V{x, a) (2.2) 



4 



N. Fernandez-Garcia and O. Rosas-Ortiz 




Figure 1. The rectangular well in a semi- harmonic background. The number of bound state energies 
Eq, . . . , E]y_i, depends on the area A — 2aVo of the well (of. Section [2.1.11) . Some of the incoming 
scattering waves of energy E > are 'captured' for a while in the region (—oo,a), as discussed in 
Section [5l If Vq~^ = 2a, the limit a — >■ leads to a semi- harmonic delta well (see Section [3]). Semi- 
harmonic rectangular and delta barriers are obtained by the change Vb — )■ — Vb (cf. Section |4]). 



is solved by the appropriate combination of the functions 

V^i(x) = Aie-"'/2 (^a, i; rr^) + B^e-'^'/^x iFi (a + i, |; x^), 
tpn{x) = A2u{x) + B2v{x) := A2 sinqx + B2 cos qx, 
V^iii(x) = A3e*'=^ + i?3e-*'", 

where k = \fE and q = y/Vo + k'^, the functions e^'^^ and e"*'^^ represent waves moving towards 
the right and the left respectively, and iFi{a,^; z) stands for the confluent hypergeometric 
function (see Appendix) with a = . 

To analyze the solutions of (j2.2p as x — )• —00, we use ()A.3p . A simple calculation gives 



iplix) 



-e 2 



(xe-*^) 



A, 



B, 



2 re(a + l/2) 



as X 



-00, 



with re(2;) the Euler gamma function of 2; G C. Therefore, to construct solutions which are 
regular at the left edge of DomF(x,a) = (—00,00), the coefficient Bi is constrained to satisfy 



B, = 2A 



re(a + l/2) 



re(a) 

As a consequence, the function ipi reads 



iFi{a,l,x^)+2x^^^^^^iFi{a + l,l;x^) 



AMx), (2.3) 



and it becomes zero as x — >■ —00. If the parabolic part of the potential appears to the right, 
rather than to the left of the well, it occurs a phase difference in the asymptotic behavior of the 
solution (compare equations ()A.2p and (jA.SP ). We decided to put the oscillator- like interaction 
to the left of the well to precisely emphasize the subtleties of the involved mathematics. 
The condition of continuity at x = —a leads to the coefficients 



A2 



(p{-a)C2, 



B2 ^ M-Pu + P^) 

Ai q 



ip{-a)D2, 



where — /3/ stands for the logarithmic derivative of the function f{x). Hence we have 



i^iiix) = Ai(p[—a)[C2 smqx + D2 cosgx]. 



(2.4) 
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In turn, from the condition of continuity at x = a, one arrives at the expressions 

Here z* is the complex conjugate of z E C, and the Jost function T{k, a) is defined as 
T{k, a) = -[^ [n(/3_ + + v{p. + P,)D2] 



k 



(2.5) 



with — /3_ the logarithmic derivative of e . Thus, for the solution in region HI we have 



^iii(x) = Anp{-a)-[F{k,a)e 



-ikx 



T*{k,a)e'^% 



(2.6) 



Scattering states. Consider now a particle of energy = A:^ > 0, arriving from +00 
towards the effective zone of the semi-harmonic well (|2.ip . The reflection amplitude 

F*{k,a) 



s{k, a) 



(2.7) 



T{k, a) 

is directly obtained from (12. Gh and has modulus one. This can be written as s{k,a) = e^^^i^'''^) , 
We see that the effect of the semi-harmonic rectangular well is to cause a phase shift (up to 
integer multiples of vr) of the reflected wave by 

lmT{k, a) 



6{k, a) 



arctan 



KeJ^{k, a) 

Then, for scattering states, the solution of (j2.2p cancels as x 
region H (see equation ()2.4p ). and behaves as 

'(piuix) = Anp{—a)\J^{k, a) I sin(A;a:: -|- 6) 

in region HI. 

2.1 Conditions for the trapping of particles 



(2.8) 

-00, oscillates harmonically in 



The solutions of the Schrodinger equation (\2.2^ which correspond to either bound or resonance 
states are picked out from (j2.3p . (j2.4p and (j2.6p such that only outgoing waves exist. In other 
words, to construct these solutions, it must be imposed the Siegert condition 



lim /3^(x) 



-ik, 



k E 



(2.9) 



This last expression is equivalent to make J-{k,a) = in (j2.6p . Then, the Jost function (j2.5p 
is analytically continued so that the reflection amplitude (j2.7p is a meromorphic function with 
poles on the zeros of J-'{k,a). For simplicity, hereafter Ai = l/ip{—a). Therefore, the wave 
functions of either bound states or resonances are written as 

ip{x) 



ip{x) 



, X < -a, 

f{-a) 

C2smqx -\- D2COsqx, |x| < a. 



such that J^{k, a) = 0. 



(2.10) 



a < X, 



-§-F*(A;,a)e^'=^, 

Given a 7^ 0, the zeros of the Jost function J-{k, a) are defined by the roots of the transcendental 
equation 



-ik 



ik(3^{-a) 



tan 2qa. 



(2.111 



It is straightforward to verify that this last equation has no solutions on the real line. 
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Table 1. The dimensionless single bound state energy of a unit area rectangular well in the free-particle 
(FP) and semi-harmonic (SH) backgrounds for different values of a. In both cases, the numerical results 
include the delta well as a limit (a — > 0). 



a 


rectangular well 


rectangular well 




(FP) 


(SH) 


2.0 


-0.113438 


-0.045272 


1.5 


-0.130400 


-0.039482 


LO 


-0.153960 


-0.033514 


0.5 


-0.189338 


-0.037435 


0.0 


-0.25 


-0.079710 



2.1.1 Bound states 



Potential (j2.ip includes a region of classical confinement in which some bound state energies 
could be present. In this context, if A; = iK is a zero of the Jost function with k > 0, then 
E = k"^ = —K? is the energy eigenvalue of a bound state. Indeed, if A; is a solution of (|2.1ip on 



the positive imaginary axis of the complex A;-plane, one has ct — 4 — 4 
The number of these roots is determined by the area A = 2aVo of the well 
bound states are represented by the square-integrable functions 



> 0, and G 



Thereby, the 



f>n{-a) ' 

C2 sin qnX + D2 cos QnX, 
u{(3u + K.n)C2 + v{(3v + K,n)D2 



X < 



-a, 



2Kr; 



|x| < a, 

g-K„(a:-a) a < X, 



n = 0,l,...,N{A) -1. 



For instance, if V^~^ = 2a, potential (j2.1|) becomes a unit area rectangular well in a semi-harmonic 
background. This admits a unique bound state, just as it is reported on Table[T]for five different 
values of a (see also Fig.[2|). For comparison, the equivalent results for a free-particle background 
are also reported. In both cases, the single bound energy goes to a definite negative number 
as a — )• 0. However, given a > 0, this energy is less negative in the semi-harmonic background 
than in the free-particle one. The same effect is observed for all the bound states belonging to 
a rectangular well of arbitrary area A ^ 0. Thus, in a semi-harmonic background, the bound 
states of the rectangular well are displaced towards the threshold. Notice that the oscillator- 
like interaction to the left of a given (symmetric) well produces a left-right asymmetry in the 
involved wave functions (see, e.g.. Fig. [2]). 



2.1.2 Resonances 

Siegert functions are solutions of the Schrodinger equation (j2.2p . associated to complex eigen- 
values e = E — i^, and fulfilling the purely outgoing condition (j2.9p . It is usual to consider e as 
a compound of the resonance position E and the inverse of the involved lifetime = F. We 
obtain these complex eigenvalues from the wave numbers k in the fourth quadrant of the complex 
fc-plane which are solutions of (j2.1ip . Some of the first resonances associated to a semi-harmonic 
rectangular well of unit area are reported on Table [2] for different values of a. In contraposition 
to their equivalents in a free-particle background, where the symmetry of the rectangular well 
is inherited to the wave functions and resonances (see e.g. [H] and [IS])) the Siegert functions 
of the rectangular well in a semi-harmonic background cancel at x = —00, obeying the left-right 
asymmetry generated by the oscillator-like interaction. This is illustrated in Fig. [3l 
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Figure 2. The wave function of the unique bound state belonging to the semi-harmonic rectangular 
well of unit area with a = 2 (continuous, blue) and a = (dashed, black). The dimensionless energies 
are Ea=2 = -0.045272 and Es = -0.079710 respectively (see Table [J). Notice the left-right asymmetry 
produced by the parabolic part of potential (|2.ip . 



Table 2. The first five resonances e = E — belonging to the unit area rectangular well in a semi- 
harmonic background for different values of a. The case a = (the semi-harmonic delta well) is such 
that Re(e) = 4m + 3 + 7m, with jm ~ 1 and m = 0, 1, . . . (see Fig. [5]). 



a 


resonances 


a 


resonances 


2 


00.623117 -iO.599545 
05.260402 - i2. 255076 
08.875649 - i2.603548 
11.217715 -il. 995880 
17.359977 - i2.229026 


0.5 


03.569260 -il.487849 
10.233701 - i2. 499541 
13.960171 - i2.280915 
17.904306 - i2. 258438 
21.878702 - i2.356575 


1.5 


01.009578 - i0.981433 
03.852457 -i2. 106492 
07.574481 - i2.017753 
09.826328 - i2. 651538 
13.723866 -i2. 196828 


5 X lO"'' 


03.792859 - i0.909297 
07.852027 - il. 117703 
11.880118 - il. 242053 
15.897157 - il.331170 
19.908829 -il.400707 


1 


01.838241 -il. 632446 
05.675163 -il. 760804 
08.749583 - i2.720863 
11.835116 -i2. 294470 
15.830298 - i2.343309 





03.792839 - i0.909196 
07.852012 - il. 117599 
11.880106 - il. 241947 
15.897146 - il. 331064 
19.908819 -il.400600 



3 The semi-harmonic delta well 



A particular case of the previous results which deserves special attention is obtained in the limit 
a —7- 0. It is well known that the family of rectangular wells of unit area converges to the delta 
well in the sense of distribution theory (see e.g. [21] and [17]). The transparency properties of 
the point interactions are such that no resonances can be associated to zizQ6{x), with Q an opac- 
ity parameter |18] . The situation changes if the free-particle background is replaced by a less 
trivial scenario. Namely, the presence of resonances is ensured for the point interactions either 
by putting a constant potential as the background [9llll] (see also [TO]), or by constraining the 
potential domain of 5{x — xq) to be [0, -|-oo) rather than the straight-line M [18]. In contradis- 
tinction, point dipole interactions represented by the derivative of the delta distribution ziz6'{x) 
admit resonances in a natural form |17] . We are going to analyze the bound states and reso- 
nances associated to a delta well in the semi-harmonic background introduced in the previous 
section. 
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Figure 3. The real (continuous, blue) and imaginary (dashed, black) parts of the Siegert function 
belonging to the first (left) and fifth (right) resonances reported on Table [2] for a = 2. In all cases, the 
left-right asymmetry is due to the parabolic part of potential ()2.1I) . 



Potential (j2.ip represents a family of functions of compact support (the unit area rectangular 
wells, parameterized by a) in a semi-harmonic background. If a — t- 0, we have 

Vs{x) := lim V{x; a) = Vs(x) — S{x), 

a— >0 

with 




a: < 0, 
x>0 



(details of the limit procedure for this kind of potentials can be consulted in [2T1 [22] ) . Applying 
the same limit to the transcendental equation ()2.1ip . one gets the expression 

2re(^)=(l + ^^)re(^^). (3.1) 

Equation ()3.1|) admits an isolated root on the positive imaginary axis of the complex /c-plane 
(see Fig. S]). We have ks = fO.2823302. The corresponding energy Es = kj = -0.0797104, 
is consistently recovered as the limit o — )• of the values reported on Table [H Interestingly, 
this bound energy is closer to zero than the bound energy of the delta well in a free-particle 
background E = —0.25. That is, the semi-harmonic background causes the displacement of the 
single bound energy towards the threshold, just as it has been remarked above. The involved 
wave function 

is depicted in Fig. [2j Observe that the parabolic part of the potential produces a decreasing 
of ip{x) which is faster than the decreasing of its counterpart e~'^^ in the free particle zone. 

In the present case, the resonances are determined by the roots of (13. ip in the fourth quadrant 
of the complex fc-plane. Up to a constant 7m, the position E^ of each of these 'complex energies' 
is clearly connected to the distribution of the odd energy eigenvalues of the harmonic oscillator: 
Re(e) = 4m -|- 3 -|- 7m, with m = 0, 1, 2, . . .. For the results reported on Table El the parame- 
ter 7m approaches 1 from below as m increases (see Fig. [5]). As a conclusion, different than the 
conventional delta well, a delta well in a semi-harmonic background admits resonances. More 
details will be given in Section [5l 



4 Semi-harmonic barriers 

The case of a particle in a semi-harmonic rectangular barrier can be analyzed as a consequence 
of the previous results. The change Vq — )■ —Vq in (j2.ip . (j2.2p produces q — >■ iy/\k'^ — Vo| for 
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Figure 4. The isolated zero of the Jost function which corresponds to the unique bound state of a semi- 
harmonic rectangular well of unit area for a = 1 (continuous, black), a = 0.1 (dashed, blue) and a = 
(continuous, red). In each case, the curve has been depicted as a function of /c = ik, with k g M. 

1 5 9 13 17 3 7 11 15 19 

OFS * i i * ^ OFS i i S 

-0.5 -0.5 

-1 - 1 -1 - 

-1.5 I 1 -1.5 I — 



Figure 5. First five resonances e ~ E — of a, delta barrier (left, red disks) and a delta well (right, blue 
disks) in a semi-harmonic background. Black disks on the real line represent the even and odd energies 
of the harmonic oscillator respectively. Vertical lines are included as a reference. The numerical values 
of e can be consulted on Tables [3] and [2] respectively. 



k'^ < Vq, and q — )■ \/k''^ — Vq for k'^ > Vq. The straightforward calculation shows that there is no 
bound states for this system. Some of the corresponding resonances, in turn, are reported on 
Table [3] (see also Fig. [6]). In this case, the complex energies e are below and close to the positive 
real axis such that, in the very limit a — ?■ 0, Re(e) is in one of the positions 4:m + l + Xm, with Am 
a constant and m = 0, 1, 2, . . .. For the results on Table [3l approaches 1 from above as m 
increases (see Fig. [5]). Then, the resonances of a delta barrier in a semi-harmonic background 
are distributed in correspondence with the even energy eigenvalues of the harmonic oscillator. 
See more details in Section [5l 



4.1 Darboux transformations 

The Siegert condition (|2.9|) is appropriate to construct the complex supersymmetric partners of 
the semi-harmonic potentials analyzed in the previous sections. Consider the Darboux transfor- 
mation 

V{x;a) = V{x;a) +2l3'{x), (4.1) 



with V{x; a) either a semi-harmonic well or barrier, and /3(x) the logarithmic derivative of any of 
the corresponding Siegert functions Condition (j2.9p indicates that the new potential V{x; a) 
behaves as V{x;a) in the limit x — )■ -|-oo. On the other hand, the straightforward calculation 
shows that (3{x) — t- 0, as x — )■ — oo. Therefore, V{x;a) is a complex potential behaving as 
V{x;a) at the edges of Doml/(x;a) = (— oo,-|-oo). The supersymmetric formalism of quantum 
mechanics [28l [23l [Ml [25l [26l [27] ensures that the energy spectrum of V{x; a) is the same as that 
ofV{x,a) ifthesuperpotential is defined as /3(x) = -£lml;,{x) (cf. [l5l[30l[3Tl[32l[33]). InFig.[3 
we show the behavior of the Darboux-deformed semi-harmonic rectangular barrier, constructed 
from ()4.ip with the Siegert state belonging to the fifth resonance of Table [3l The presence of 
maxima and minima, in both the real and imaginary parts of the potential, makes V{x; a) to 
behave as an optical device which both refracts and absorbs light waves [231 [IS] • A similar result 
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Table 3. The first five resonances e = E — belonging to the unit area rectangular barrier in a 
semi-harmonic background for different values of a. The case a = (the semi-harmonic delta barrier) is 
such that Rc(e) = Am + 1 + A,„, with A,„ ^ 1 and m = 0, 1, 2, . . . (see Fig. [5]). 



a 


resonances 


a 


resonances 


2 


00.595222 - iO. 312336 
04.211813 -il. 853317 
07.387244 - i2. 485184 
10.572524 -il.985917 
12.865432 - i2. 874092 


0.5 


01.957470 - i0.802720 
05.720407 - il. 561124 
09.297545 - i2.218137 
12.626496 - i2. 505935 
16.213699 - i2. 433756 


1.5 


00.846017 - iO. 497624 
06.188284 -i2. 213815 
08.833155 -i2.011628 
12.476467 - i2. 723854 
14.938732 - i2. 618887 


5 X 10""^ 


02.076264 - iO. 718026 
06.065577 - il.025854 
10.058023 - il. 181640 
14.052696 - il. 286283 
18.048705 - il.360061 


1 


01.313699 - i0.776043 
04.317689 - il.833635 
07.338936 - i2.031873 
10.901777 - i2. 284989 
14.229903 - i2. 730174 





02.076211 - iO. 718123 
06.065549 -il.025956 
10.058003 - il. 181793 
14.052679 - il. 286389 
18.048690 - il.365109 




Figure 6. The real (continuous, blue) and imaginary (dashed, black) parts of the Siegert function 
belonging to the first (left) and fifth (right) resonances reported on Table [3] for a = 2. Remark that the 
left-right asymmetry of all these functions is generated by the parabolic part of potential (|2.ip . 

is obtained for the Darboux-deformations of the semi-harmonic wells. The procedure can be 
repeated at will, giving rise to more elaborate complex deformations of the initial potentials. 
Particularly interesting, the double transformation obtained by using a complex eigenvalue e in 
a first step, and its complex conjugate e* in the second step, produces deformations which are 
real functions (see, e.g. |3H I15 1 [7]). In such a case, the semi-harmonic barriers will exhibit 'hair' 
as a characteristic of the deformation [15^ [7] . 

5 Dwell time 

As yet the bound states and resonances have been analyzed as different manifestations of the 
same sort of mathematical solution. Both of them are defined in terms of the zeros of the Jost 
function (j2.5p and are represented by a Siegert state of the form (j2.10p . Concerning scattering 
states, let E > be the energy of a particle which impinges on a scatterer from the right, 
according to the rule (j2.ip . Here Vq can be either positive or negative. Exterior to the inter- 
action zone (— oo,a), the stationary scattering wave function is -0111 g"*'^^ — s{k,a)e^^^ (see 
equations p.6|) and ()2.7p ). Since there is neither sources nor sinks we have total probability 
conservation |s(A;,o)p = 1. We want to get some insights on the time spent by the incident 
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Figure 7. Left. The real (continuous, blue) and imaginary (dashed, black) parts of the supersymmetric 
partner V{x; a) of a semi-harmonic rectangular barrier. This has been constructed with the fifth resonance 
for a ~ 2. Right. The Argand diagram of V{x; 2). The circle o and disk • serve as reference and correspond 
to V{—3; 2) and V{3; 2) respectively. 

particle in the interaction zone. To analyze the 'capturing' of the particle by the semi-harmonic 
rectangular potential, we assume |s(A;,a)|2 = e-^S/-/, with rj?, the mean time spent by the par- 
ticle in the interaction zone (—00,0), and Tf a characteristic time constant to be determined. 
Using the absorption probabilities method \50\ [5T| 152] , a small imaginary part AVj is added 
to the potential under consideration. The absorption dimensionless Schrodinger equation to be 
solved is {H + iAVi)ipa = i^- 

Given a solution ipo of the unperturbed equation (for which 
AVj = 0), we use the ansatz V'a = ipo9 to get V'a = ipoe^^'* (an arbitrary integration constant 
has been omitted for simplicity). Therefore, the probability density of finding the particle is 
l^ope^^^^*. The time Tf required for the probability density to decrease a factor of its initial 
value IV'oP is Tf = — 1/(2AV7), as this is introduced in [51]. The dwell time (or sojourn time) 
Tf) in the region (—00, a) can be now obtained from |s(A:,a)p = e^'^S^^^. Indeed, in the limit of 
small AVj, a simple calculation gives 

T?j = - lim -^|s(A;,a)|2. (5.1) 

Now, consider a solution of the absorption eigenvalue equation {H + iAVi — e)il)e = 0, with 
e £ C. The derivative with respect to AVi of this last equation, after evaluating it for AVi = 0, 
is reduced to 



dAVi 



= 0. 

AV/=0 



The derivative of the absorption eigenvalue equation with respect to e/, the imaginary part of e, 
after evaluating it for AVi = 0, produces 

dei 

Adding these last two results yields 

a d 

lim — = — — — . 

AVf^o dAVi dei 

In this way, equation (j5.ip can be rewritten as follows 

Th = -l^^\s{k,a)\'. (5.2) 

Since = e, the reflection amplitude s{k, a) must be a complex analytic function of the complex 
eigenvalue e = E + iej. The Cauchy-Riemann condition reads as 

dsB__ds]_ dsR _ dsi 

dE ~ dei dei ~ dE " 
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The introduction of these last expressions into (|5.2p produces 

= 2^ = rw, (5.3) 

where we have used s{k, a) = e^^^^^''^^ . From the above expression, the dweh time of the scattering 
particles which are 'trapped' in the semi-harmonic rectangular potentials is in direct connection 
with the slope of the reflected wave phase shift S{E, a). Thereby, the dwell time ()5.3p coincides 
with the definition of phase time Tyy, introduced by Eisenbud and Wigner as the time delay in 
binary collisions [36]. In particular, the formation of a resonance e = E^ — i^ introduces a positive 
time delay between the arrival and departure of the scattering particle from the region (— oo, a). 
This last means that particles impinging the semi-harmonic rectangular potentials with energy 
E = Er, will spend times in (— oo,a) which are larger than the times spent by the scattering 
particles of energy E = Er + e, with e ^ a real number small enough. A rapid increasing of 
the phase shift is then expected in the vicinity of the resonance position Er- Hence, the phase 
shift (I2.8p encodes enough information to identify true resonances by associating the peaks 
of ()5.3p with the real part of the complex eigenvalues e (see, e.g. [M]). Of special interest, the 
time delay of point interactions in a semi-harmonic background can be calculated as the limit 
o — >■ of (15.31) . One obtains 



- \\^Wl[E) ) -^E- ^^-^^ 
where the sign + (— ) stands for the semi-harmonic delta barrier (well) and 
Im(J"±) -Te (^) \/^ 

An straightforward calculation shows that — )• ^, as S — +oo. From Fig. [HI we notice that 

Ty^^ (r^^) exhibits a series of Fock-Breit-Wigner peaks Wl, W2, ... (Bl, B2, ...), centered 
each one at the real part of the first, second, etc, energy resonances reported on Table [2](Tab- 
le[3|). Thus, the zeros of the Jost function (j2.5p are interrelated with the maxima of the time 
delay (|5.4p . as expected for true resonances [3l]. A conclusion which makes self-consistent our 
approach since the same results are obtained in two different forms. From Fig. [HI we also notice 
that all the 'captured' particles spend a time tw > 1 in (— oo,a), except for a region of small 
positive energies in the rectangular well case. The closer to the resonance position Re(e) is the 
incoming energy E, the larger is the time spent by the particle in the interaction zone. 

In general, after entering the region (—00,0) of an arbitrary semi-harmonic rectangular po- 
tential, the particle may be reflected back and forth many times before it manages to get back 
out (see discussions on the matter in [35]). These reflections may take place either in the zone 
of the rectangular potential (—a, a), the parabolic part of the potential (—00, —a), or in a com- 
bination of both. Our calculations consider the complete scattering process so that it is not 
possible to distinguish among all possible internal reflections and self-interference during the 
approach to the region (—00, a). For an arbitrary isotropic potential Vq (i.e., Vq{x) becomes 
zero as |x| — )■ +00), the dwell time td must obey the identity td = trI-RP + rr|Tp, with R{T) 
the reflection (transmission) amplitude, and tr {tt) the average time spent in Vq by reflected 
(transmitted) particles. For instance, if Vq is a rectangular barrier in the free-particle back- 
ground, one has tr = tt = td, so that all particles spend, on average, the same amount of time 
in Vq, no matter the size of the barrier [531 [13 [50] . Our expression ()5.2p is then associated to the 
reflection time, the mean time if particles are finally reflected, in the absorption model |51j (see 
also [53]). The comparison of r^, as this has been defined in (|5.3p . with tr for the rectangular 
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Figure 8. The Eiscnbud-Wigner time delay tw for the scattering process of a semi-harmonic delta 
barrier (dashed, black) and well (continuous, blue). The peaks Wn (Bn) localize the real part of the 
resonances reported on Table [5] (Table [3]). 

barrier could give us additional information of the time delay induced by the semi-harmonic 
background. 

Roughly speaking, time delay of the semi-harmonic rectangular potentials should correspond 
to the excess time that the scattering particle spends in (— oo, a), when compared to a free particle 
under the same initial conditions. Since our model considers the semi-harmonic background as 
the environment into which the square potential is embedded, this last can be taken as an open 
quantum system [5]. In this respect, the above indicated comparison of and tr makes sense 
in order to calculate time delays as the difference of reflection times. Some insights would be 
also obtained concerning transmission times. Given an open interval lo = {xi,X2) of sufficiently 
large size X2 — xi in (— oo,a), and a time window 6 = (ti,t2) with t2 > ti, one could introduce 
the corresponding dwell time as 



with H as this has been given in ()2.2p . and an scattering state of energy E. In a similar form, 
the expression 



with Hq the Hamiltonian of the square potential in a free-particle background, could be taken 
as the dwell time of the 'free-particle' system. The time delay in the space-temporal window 
(oj, 9) should correspond to the difference — Go, as usual. In such a definition the quantities Q 
and 9o are assumed to be well defined and finite for arbitrary uj G (— oo,a). In principle, it 
seems to be the case since all the involved scattering wave functions cancel at x = — oo (this 
is because the parabolic part of the environment). In order to arrive at a quantity which is 
independent of w, one could take the limits xi — )• — oo and X2 — ?• a. However, the anisotropy of 
the semi-harmonic rectangular potentials makes this last step not evident a priori. That is, some 
caution is required to apply the difference of sojourn times as the definition of time delay in 
the semi-harmonic rectangular potentials. Mainly, as connected with transmission dwell times 
in the interaction zone (— oo, a) since the potential diverges as x'^ at x = — oo. Thus, in this last 
case, there is not a clear connection between sojourn times and phase times because almost all 
the approaches on the matter require potentials localized in a finite zone of the real line. The 
problem deserves a detailed analysis and the results will be reported elsewhere. 
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6 Concluding remarks 

In this paper, we studied the energy properties of a particle in the presence of a semi-harmonic 
rectangular potential. The latter is a one-dimensional rectangular potential in a background 
composed by an harmonic interaction to the left and a free particle interaction to the right 
of the rectangular potential. The advantage of this model lies in its simplicity to get analytic 
expressions for the wave functions, scattering states and resonances. Indeed, the problem is 
faced by considering bound states and resonances as different manifestations of the same sort of 
mathematical solution. Both of them are defined in terms of the zeros of the Jost function and 
are represented by a Siegert state. The numerical integration of the energies is done by solving a 
simple transcendental equation. In particular, it is found that, in the presence of a semi-harmonic 
background, the conventional bound energies of the rectangular well are displaced towards the 
threshold. Such behavior is preserved in the limit where the well becomes a delta. In this limit, 
there is a single bound state of energy Eg = —0.079710, which is less negative than its counterpart 
in a free-particle background, the value of which is E = —0.25. The resonances appearing for 
the semi-harmonic delta well exhibit a very peculiar behavior. They are distributed below and 
close to the positive real axis of the complex plane in such a way that their real part mimics the 
odd energy eigenvalues distribution of a harmonic oscillator. Namely, they are located according 
to the rule 4m -|- 3 + 7m, with 7.^ ^ 1 and m = 0, 1, 2, . . .. A similar situation occurs for the 
resonances of a semi-harmonic delta barrier, with the distribution ruled by 4m + 1 -|- A„ with 
Xm ~ 1- Thus, in this case, the semi-harmonic background induces the real part of the resonances 
to be distributed in correspondence to the energy eigenvalues of the harmonic oscillator. The 
Siegert functions were successfully applied to construct complex Darboux-deformations of the 
semi-harmonic rectangular potentials. These new potentials behave as optical devices which 
both refract and absorb light waves. 

Concerning the times involved in the scattering process, we assumed the time exponential 
decay rule \s{k,a)\'^ = exp(— r^/rj), with s{k,a) the reflection amplitude, rj a constant to be 
determined and the mean time spent by an scattering particle coming from +00 towards the 
interaction zone (—00, a). The application of the absorption probability method [501 151^ [52] al- 
lowed us to show that coincides with the phase-time of Eisenbud and Wigner tw, calculated as 
the derivative of the reflected wave phase shift 6{E, a) with respect to the energy E. Becuase the 
global properties of the potentials we have discussed on, the dwell time so calculated corresponds 
to a reflection time (the mean time if particles are finally reflected) . In the limit where the rect- 
angular potentials become point-like interactions, we derived an explicit analytical expression for 
Tw, and showed that its local maxima are in correspondence with the resonances obtained in the 
previous sections. Finally, the semi-harmonic rectangular potentials can be viewed as the open 
quantum systems integrated by a rectangular potential (the system itself) and an environment 
(the semi-harmonic background). In this picture, the system is an isotropic potential for which 
the time spent in any finite region of space, averaged over all incoming particles, is well defined 
and corresponds to the dwell time td involved. Therefore, the semi-harmonic rectangular poten- 
tials can be seen as the 'system' affected by a semi-harmonic interaction. Some insights on time 
delay can be obtained by calculating differences between the dwell time of these two systems 
in a given space-temporal window. This approach, however, is not directly applicable in the 
present case because the anisotropy of the environment. Work in this direction is in progress. 

A The confluent hypergeometric operator 

The kernel /C(a,c) of the confluent hypergeometric operator 




— a. 



a,c G M, 



z G C 
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Figure 9. The two first sheets of the Riemann surface defined by the factor e in equation (|A.1|) . the 
cut is at -1. If < argz < then and 6 = e~^'^*z^ are in the second and first sheets respectively. 

is integrated by the solutions of the Kummer equation L(^a,c)f{0'^ c; z) = 0. That is, if /(a, c; z) is 
a confluent hypergeometric function then f{a,c;z) G K,{a,c) [MIES]. The asymptotic expansion 
of f{a, c; z) for large values of \z\ is given by 

^) = f^^^^""^~"[l + + f|^"'^""'[^ + - «)] (A.l) 

where the positive sign in the factor e^"'^* is taken when — f < argz < while the negative 
sign is taken when — ^ < arg z < ^ [56] and 

P±(7) = f;(±r ^^^"^^J^'^^^" , (7)n = 7(7 + l)(7 + 2)---(7 + n-l), (7)0 = 1. 

n=l 

The above expressions can be utilized to deduce the asymptotic expansion of f{a, c; z) for specific 
ranges of argz. For instance, if — ^ < argz < ^ we get 

/(a,c;z)^^e^z"-^[l + P+(c-a)], Re(z) > 
r(a) 

and 

f(a,c;z^)^l^e'\^^''-'^ if z ^ +00. (A.2) 
r(a) 

To find the behavior of f{a,c;z^) for z — )• —00 we consider a range of argz which includes 
the negative real values of z. For this, let us assume ^ < argz < ^ (see Fig. [9]). Then 
^ < argz^ < and z^ is in the second sheet of the Riemann surface defined by the factor 
g-avri (jA.ip . Since = e'^^'^z^ is such that < arg6' < |, from (jA.ip we arrive at 

f{a,c;9) = ^^^""^ . e°"z-^"[l + P_(a)] + lMe^^2('^-^)e-^"(°-^) [1 + P+(c - a)]. 
L [c — a) 1 (aj 

Therefore, if arg z = vr one finally gets 

/(a,c;z2) I|^e^'z2(a-c)g-2^i(a-c) if z^-cx). (A.3) 
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